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We give predictions for the neutrinoless double beta decay rate in a simple variant of the A4 
family symmetry model. We show that there is a lower bound for the PPo^ amplitude even in the 
case of normal hierarchical neutrino masses, corresponding to an effective mass parameter |mee| > 
0.17\/ AmATM- This result holds both for the CP conserving and CP violating cases. In the latter 
case we show explicitly that the lower bound on \m^e\ is sensitive to the value of the Majorana 
phase. We conclude therefore that in our scheme, PPov may be accessible to the next generation of 
high sensitivity experiments. 

PACS numbers: ll.30.Hv, 14.60.Pq, 12.60.Fr, 12.60.-i, 23.40.-s 



Are neutrinos their own antiparticles? This is one of 
the most basic current unknowns in neutrino physics. Of 
all possible manifestations of neutrino masses, neutrino- 
less double beta decay [PPqi, , for short) offers - to date 
- the only potentially viable way to answer this question. 
If /J/Sof exists, then neutrino masses are Majorana in na- 
ture, irrespective of their ultimate origin m. Together 
with cosmology 0] and direct kinematical searches in tri- 
tium decay [3, PPqi/ offers one of the three main com- 
plementary ways to probe the absolute scale of neutrino 
masses. 

Current experimental limits from /?/3oiy on the effec- 
tive Majorana mass of the neutrino m^e are of order 
rUee < 0.3—1 eV 0, We note that a claim for a 
finite /3/9o!y rate has been published in j^, but this has 
so far not been confirmed by any other experiment. A 
recent proposal [?,] aims explicitly at testing the half- 
life range suggested in (gj. However, future /3/3oi/ exper- 
iments may be able to reach down to much lower mass 
scales. Many experiments sensitive to — 0.05 eV 
have already been discussed; see for examp le la, Ip . In 
the longer-term future, even rriee — 0.01 eV 'io^s 
not seem impossible. 

The historic confirmation of neutrino oscillations over 
the last few years [l^, together with some basic the- 
ory, suggests that /3/3oi/ is expected, although in general 
no lower bound on the magnitude of the expected ef- 
fect can be given. Theoretical input is therefore needed. 
Currently the origin of neutrino masses is completely un- 
known. The basic dimension-five operator which leads 
to neutrino masses jl3j can arise from a variety of mech- 
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anisms characterized by vastly different scales. Alter- 
natives include the seesaw mechanism lisl I16 I and 
low-energy i?— parity violating supersymmetry [ij]. In 
neither case is it possible, in general, to establish a lower 
bound on the magnitude of the /9/3oi/ rate. 

Here we consider a simple phenomenological model 
based on a new realization of the family symme- 
try El El H in which no SU(3) ® SU(2) ® U(l) sin- 
glet neutrinos are introduced. Instead, the small neutrino 
masses arise from the small induced vacuum expectation 
values (VEVs) generated for the neutral components of 
triplet Higgs bosons [isL W^ . transforming nontrivially 
under the family symmetry. The lepton and Higgs 
particle content and their transformation properties un- 
der Ai and SU(2) ® U(l) are specified in Table I. With 
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TABLE I: Lepton and scalar boson quantum numbers 

these transformation properties, the charged lepton mass 
matrix is already diagonal in the flavor basis, with 

TOe — hiVi + h2V2 + h^Vs 

rrin = hiVi + ujh2V2 + uj^ h^vs 



hlVl + U! h2V2 + Ojh3V3 



where hi are charged lepton Yukawa couplings, Vi — ((/)^) 
and a; is a complex cubic root of unity satisfying 1 + a; + 
uj^ = 0. The neutrino mass matrix is then of the form 

a + b + c f e 

Af„ = I / a + ujb + uj^c d I (1) 

e d a + ui'^b + loc 
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where the only non-diagonal entries are those of the A4 
triplet ^. Here we have defined 



Ai(^?> 



a 

b — ^2\'l2l 



(2) 



where Ai, k and (f?°), 
and VEVs, respectively, 
the conditions h = c, and 
the former is an ad hoc 
be maintained naturally 
symmetry. Note that M^, 
possibly unique) property 
renormalized by the charj 
in the same way, i.e. \hi 
being proportional to + 



are triplet Yukawa couplings, 

Let us further assume that 
d = e — f hold. Whereas 
assumption, the latter can 
because of a residual 
has a very remarkable (and 
here in that each entry is 

;ed-lepton Yukawa couplings 
^ + |/i2p + |^3p, instead of 

71? as in the Standard Model. 



It is easy to see that in the above limit we have the 
prediction 



723 



7r/4, 







which matches well with the neutrino oscillation 
data 12J. Furthermore it can be shown that, in the limit 
where the solar mass splitting is neglected, b and d can be 
made real, so that the atmospheric neutrino mass split- 
ting takes on a very simple form 



Am32 ~ 6bd = Am^,^ 



(3) 



The solar neutrino mass splitting Amg 
be expressed as 



SOL < ^wItm can 



Am^i = jTf + Ti + T.^ = Amt, 



where 



Ti = 672151 |d|sin((/.2) 

T2 = 2V2|d|(2|a|cos(</>i 

T3 = -3|6p + - 6|a||6| cos(0i -t- (j)2) 
+ 2\a\\d\ cos((/)i) - 2|6||d| cos(02) 



(4) 



(5) 



|5| COS(02) + Ml ) (6) 
(7) 



with <pi = (pa ~ (pd and (l>2 = 4>d — 4>b, where a — 
|a| exp(z(/<a), b — |5| exp(i(/)b) and d = \d\ expiicpd) . The 
condition in Eq. Q leads to three inequalities 

m\<Aml, 

which, normalized by Ar7i^.j,^j, can be expressed in terms 
of the small parameter a = ArngQ^/^'^iATM 



V2 



V2|sin(02)| < 



2|a|cos((/)i) -I- \b\cos{(t>2) + \d\ 



< a 



(8) 



(9) 



m\d\ 



3\bf + \df ~ 6\a\\b\ cos(?ii + 02) (10) 



-h2|a||d|cos((/)i) - 2|&||d|cos(02) 



< a 



where the current allowed values of a are shown in 
Fig. 13 of Ref. [13. 



As for the solar mixing angle, it is given here by 

2^2^ 



t2S = tan(20i2) = 



3b -d 



which reduces to 

in two ways, namely, 
6 = 0, 



tan^ 6*12 = 1/2 



b = 2d/3 



(11) 



(12) 



(13) 



Current fits of solar, reactor, atmospheric and accelerator 
neutrino oscillation data lead to a best fit point to the 
solar mixing angle for which tan^ 6*12 is slightly less than 
1 /2. Performing a series expansion of tan^ 9i2 around the 
two solutions in Eq. l(T^ . we get 



tan^ di2 
tan^ 9i2 
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(14) 
(15) 



These two branches are depicted in Fig. ^ The positive 
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FIG. 1: Solar mixing angle, tan'^^i2, vs. hjd. 

horizontal axis corresponds to a normal hierarchy (NH) 
neutrino spectrum, while the negative part corresponds 
to inverse hierarchy (IH), as depicted in Fig. This 
behavior is also recognized in the model of |i2Z |. 

We now turn to neutrinoless double beta decay. The 
neutrino-exchange amplitude for this process is simply 
given by 

{vTLy) — niee = a + 2b (16) 

In the case of real parameters, we can solve the following 
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system of equations 



= Ato^^ = |2a + 6 + 36)2 + 8^2 
= 6bd 
2V2d 



Am32 



t2s = tan(20i2) 



36- d 



and express the parameters a, b and d in terms of ex- 
perimentally measurable ones Amgo^j Arn^^^ and ^25- 
Substituting in Eq. we can therefore express {rrii,) 
in terms of these measured observables. We then obtain, 
up to an overall sign, 



Sign [Am 
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atmJ 
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(17) 



2S 



±Sign[2\/2i2S + 1 



2S 



''2S 



2S\ 



The calculated values of |mee/ yAmJ^j as functions of 
t2s according to Eq. p7|) are shown in Fig. [3 It can 




FIG. 2: |mee/V'^"^ATM| vs. \t2s\ according to Eq. (QTl. The 
bound depends slightly on the value of a: the left panel cor- 
responds to Of = 0.022 and the right one to a = 0.065. The 
dark (red) lines correspond to normal hierarchy, while the grey 
(green) line is for the inverse hierarchy case. The vertical line 
corresponds to the horizontal line in Fig.^at tan'^ 9i2 — 1/2. 

be seen that given the currently allowed experimental 3cr 
range 



tan2 6ii2 e [0.30, 0.61] 



we can set lower bounds for |TOee/-\/AmJp^|: 



0.17 < \m,J^Aml^J for NH 

for IH 



0.70 < Ime 



(18) 



(19) 
(20) 



Note that the lower bound on /3/3oiy for the NH case is 
especially relevant, as it is absent in the generic case. It 
results from this specific realization of the A4 symmetry 
and is also in contrast with previous A4-based models 
that led to quasi-degenerate neutrinos 0, [^^ ■ 

Moreover, should future precision experiments narrow 
down the experimental range for tan^ 0x2 then we might 



be able to distinguish between both neutrino mass hi- 
erarchies. For example, if tan^ 6*12 < 1/2 could ever be 
established, then we would have 



0.23 < \mee/\/Aml^J < 0.41 for NH (21) 
0.70 < \mee/y/Aml^J for IH (22) 

It can also be seen that, up to order a corrections, rriee 
can only be zero if tan^ 0x2 = 1, now strongly rejected 
experimentally. Note that the two solutions in Fig. |21 
correspond to the two branches depicted in Fig. ^ One 
can see that only in the branch corresponding to the so- 
lution 6 = and for values of tan^ 012 that are less than 
1/2, there is a relative minus sign between 6 and d, which 
is the condition for inverse hierarchy, as can be seen from 
Eq. ®. 

In the general case of complex parameters, lower 
bounds on |TOee/\/Am2^| can also be established for 
each hierarchy. This task is simplified greatly by taking 
into account the reliable approximations sm(f)2 — 0, and 
2|a|cos(/)i -I- |6|cos02 + Ml —0. From Fig. we see that in 
the complex case, lower bounds on |?nee/-\/Am2^^,| are 
indeed also established for each hierarchy. By compar- 
ing Fig. O with the right panel in Fig. [21 we find that 
these lower bounds are in fact exactly the same as ob- 
tained in the real case. The robustness of these bounds 
is easily understood. It follows from the fact that the 
maximum degree of destructive interference between the 
three neutrino-exchange contributions occurs in the real 
case with appropriate CP parities [23 |. 
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FIG. 3: Imee/^A m^TMl vs. 1*25 1 for a = 0.065 when complex 
parameters are allowed. 

The left-rising lines correspond to inverse hierarchy, 
while all others refer to the case of normal hierarchy. The 
lower-lying pairs of solid lines correspond to the lower 
bounds already discussed in Fig. El Again, the vertical 
line corresponds to tan^ 0i2 = 1/2. Finally, the remain- 
ing lines correspond to non-zero values of the relevant 
CP-violating phase varying cos((/)i) over the range [0,1] 
in equally-spaced steps. It is conceptually interesting to 
note that this phase is "Majorana type" 0,113, as we 
are still considering the case 613 = 0. 
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Last, but not least, we find it very interesting also 
that the lower bound on |mee/ y^Arnf^j which we have 
obtained depends on the value of the Majorana violat- 
ing phase |cos(0i)|. We can see from Fig. 0] that the 
lower bounds on \meel \/ Ato^^^, | for each hierarchy be- 
come weaker for cos(0i) = 1. Finally, we mention that. 




'•\) 0.2 0.4 0.6 0.8 1 
I cos((|)j) I 



FIG. 4: Lower bound on |mee/\/ lS.m\^^\ vs. | cos(<j!>i)|, given 
current t^s uncertainty. The lines in dark (red) and grey 
(green) correspond to normal and inverse hierarchy respec- 
tively. 



even though for simplicity we have focused on the case 
b — c, all results can be generalized to the case & ^ c, in 
which case ^13 is allowed to be nonzero. 

In summary, we have given predictions for the neu- 
trinoless double beta decay rate in a simple hierarchical 
variant of the A4 family symmetry model. We showed 
that there is a lower bound for the /3/3oi/ amplitude even 
in the case of normal hierarchical neutrino masses. We 
have seen that the bound is robust as it holds irrespective 
of whether CP is conserved or not. In the latter case we 
show explicitly how the lower bound on |mee| is sensitive 
to the value of the Majorana phase. Our scheme suggests 
that neutrinoless double beta decay may be within reach 
of the next generation of high sensitivity experiments. 
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